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PREFACE 
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the United States Naval Postgraduate School, Annapolis, Maryland, 
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Science in Aerology, The idea of the investigation was suggested by 
Dr, F, L. Martin of the staff of the Department of Aerological 
Engineering of the United States Naval Postgraduate School, His 
assistance and guidance in the preparation of this thesis is grate¬ 
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CHAPTER I 


INTRODUCTION 

The problem of the movement of surface lows and the forecasting 
of such movement has long been a paramount problem to the forecasting 
meteorologist. To date no infallible technique has been developed to 
accomplish this aim. No successful theoretical or statistical formula 
has been set forth so that this complex physical process may be objec¬ 
tively predicted. Success in most instances is the result of long 
experience and intelligent and thorough analysis rather than from mas¬ 
tery of "weather 11 by the tools of mathematics and physics. 

The complexities of the processes involved in measuring meteor¬ 
ological phenomena are well known. However investigations continue 
in an attempt to solve the problem, or to shed whatever light possible. 
The purpose of this paper is to investigate the relationship of momen¬ 
tum changes at the 500 millibar level to the accelerations and decel¬ 
erations of surface low pressure centers. The data for this investi¬ 
gation has been made possible by the publication recently of Northern 
Hemisphere synoptic 500 millibar and sea level charts (5). 

The theory of the general circulation of the earths atmosphere 
as set forth by Rossby (3) states that surface lows are tertiary 
phenomena imbedded in the west to east flow of the zonal westerlies. 
Such cyclones are the ones considered in this paper. Rossby points 
out that the westerlies in the middle latitudes form part of a reverse 
meridional circulation cell which is driven by the direct cells to the 
north and south. If the surface of the earth were homogeneous and the 
energy received from the sun constant, the circulation of the earth’s 



atmosphere would be simple and we would have no turbulent eddies in 
the zonal westerlies. However the earth’s surface is not homogeneous 
and the incoming solar radiation is not constant. These two factors, 
and possibly others (for example, friction) cause a break-down in the 
stability of the general circulation, and as a result the pattern of 
migratory cyclones, anticyclones, troughs, and wedges observed on the 
daily weather chart develops. As Haurwitz states (4) these turbulent 
eddies in the westerlies cause a meridional transport of heat and mois¬ 
ture . They also cause precipitation and outbreaks of cold air and so 
are of prime importance. 

These outbreaks of cold and warm air at the upper levels may be 
seen in relation to the type of circulation taking place. At times 
the wind blows strongly from west to east with little north-south 
component. At such times outbreaks of cold and warm air are diminished. 
At other times there are pronounced wedges and troughs on the upper air 
chart. Then the north-south components of the flow are strong and the 
outbreaks of cold and warm air is pronounced. 

The forecasting of the movement of these troughs and wedges is 
necessarily important to the forecaster. A formula for this purpose 
has been developed by Rossby (1 ). 

It is 

PL~ 

(1) C = 


where C = speed of trough 

U - mean west wind component 
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p m rate of change of Cariolis parometer with latitude 
L = wave length of trough 

r 

This formula has been used with some success (1) by the U. S. Weather 
Bureau long range unit. 

A period of strong west-east flow is known as high index. A rela¬ 
tively weak west-east flow is known as low index, in such a case it is 
generally considered that the north-south flow will be relatively strong. 
However this definition is not comprehensive enough to adequately de¬ 
scribe all flow patterns. For example let us take a situation with 
weak west-east flow and little or no north-south component. This is a 
low index situation but very different from the normal concept of low 

index. Another case is that in which there is very strong northwest- 

# 

southeast and southwest-northeast flow. This by index measurements is 
high index, but it has a large north-south component. 

A new, more complete idea of index has been put forth by Willett 
(6). He not only considers east-west (zonal) but also north-south 
(meridional) index (6). With the use of these two measures of index 
the flow pattern can better be objectively described. This investigation 
uses the concepts of both zonal and meridional index, and attempts to 
compare the daily changes of each with the daily accelerations of sur¬ 
face lows. 

Actually the day to day changes in zonal index are generally used 
subjectively as a guide to accelerations and or decelerations of pressure 
systems in view of (1). However since zonal index is calculated from the 
belt 35-55° N. latitude all around the globe, it could be, and very often 
is, a relatively coarse index when applied to axy individual migratory 
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low. Thus the basic idea of this thesis is to develop more sensitive 
indices than either zonal or meridional, but analogous to these. To 
this end it was decided to examine the 24-hour changes of specific 
momentum at the 500 millibar level immediately above the sea level low 
and at systematically determined points nearby in the 500 millibar 
windstream. The relationship of these points to the sea level low 
will be shown in Figure 1. 

Certain similarities in concept may be noted between this approach 
and that of ”steering”. Numerous investigators including Austin (7), 
Longley (8), and Runk (9) have found varying degrees of correlation of 
the movement of the sea level low with the actual wind velocities at 
several upper levels including the 500 millibar level. However, it is 
to be empphasized that the present approach may be thought of as an 
attempted refinement of such investigations inasmuch as accelerations 
of the sea level low are to be correlated with changes in specific 
momentum. It is to be noted that both of the last two variables are 
essentially differences of variables discussed by the above mentioned 
authors. 
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CHAPTER II 


THEORETICAL CONSIDERATIONS 

On the 500 millibar chart we will consider this distribution of 
specific momentum in the vicinity of the sea level low. Specific 
momentum is defined as pv where v is the vector velocity of the 
wind field at any point, and p> is the density at this point. From 
the equation of state, we have 


= 

where p = pressure in millibars 
T - temperature in °Kelvin 
Ra = gas constant per unit gram of air. 

Since we are dealing with a constant pressure chart we find that 

, x-i: 

p 

the value of H - ^ being a constant for the 500 millibar level. 
The value of the specific momentum is then effectively determined, 
apart from the constant K, by the simple relation 





Five values of the specific moment vim on the 500 millib ar chart 
are considered for every position of a sea level low as indicated in 
Figure 1. The five values are taken at south, east, north, west, and 
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We will investigate the effect of changes of specific momentum in 
these sections taken singly and in various combinations upon the accel¬ 
eration of the surface low pressure system. 

The data was taken from the historical weather map series published 
by the Army Air Forces (5) using the 500 mi Hi,bar and sea level charts 
for the same day. Unfortunately there is a time discrepancy of 8g hours 
between the map times of the two charts, the map time of the surface 
chart being 1230 GMT and that of the 500 millibar chart 0400 GMT. 
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Observations were taken only in the Northern Hemisphere, but otherwise 
there were no restrictions as to the location of the surface lows 
considered if 500 millibar data were available. Each individual low 
was followed until a central position could no longer be defined on the 
map. The observations were taken for the months of October and November 
1945. 

The unit of velocity used is degrees latitude per day (°Lat./day) 
which is numerically equal to miles per hour multiplied by 0.338. 
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CHAPTER III 


STATISTICAL RELATIONSHIPS DEDUCED 
FROM THE DATA 

The methods of collecting the data and computing the statistical 
variates •w ill be described in some detail to lend clarity to the method 
of attacking the problem and the meaning of the results obtained. 

In general lows that maintained their identity for 4-5 days or 
longer were chosen for observation. When a surface low was first ob¬ 
served on the chart, its position and the specific momentum values at 
the 5 sections (see Figure l) were recorded. One day later the position 
of the low was again recorded and its northward and eastward components 
of motion noted. The 500 millibar specific momentum components were 
recorded section by section and for each section the 24 hour differences 
in these components were recorded. These differences in specific 
momentum components were also recorded for each subsequent day that the 
low was observed. 

On the third day, two successive displacements of the surface 
center have been recorded. By taking the difference in these displace¬ 
ments (eastward and northward components) we obtain essentially the 
eastward and northward acceleration components of the surface low. 
Acceleration components are also computed on all subsequent days until 
the position of the low is no longer clearly defined. 

It is readily seen that for the calculated acceleration components 
which exist on any particular day, there are available values of specific 
momentum changes both for that day and the day before. This enables us 
not only to compare concurrent specific momentum changes and acceler- 
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ations of the surface low, but also to compare the acceleration of 
the surface low to the specific momentum changes which occurred 1 
day previously. The specific momentum changes concurrent with the 
acceleration of the surface low are referred to as "no lag" or "not 
used for lag". The specific momentum changes which precede the 
surface low accelerations are referred to as "used for lag" or "lagged 
1 day". 

The temperatures observed ranged from 268° K. to 228° K. or 40° K. 
Thus the values of specific momentum could differ by as much as 17 • 5% 
from values obtained if temperatures were not considered. Thus it is 
felt that the difference in temperature is a significant factor. 

As stated in Chapter II all velocities are converted into °Lat./ 

day which is a convenient measure of the velocities of the surface low. 

Since the magnitudes of the differences in specific momentum are of the 
-2 v 

order of 10 in the units of t employed, their values are multiplied 
by 100 to make them of the same order of magnitude as the accelerations 
of the surface low. 

Such data were recorded and computed for a number of low pressure 
systems appearing on the October and November 1945 weather maps. The 
arithmetic means and standard diviations were calculated from the 
standard formulas shown below and are tabulated in Tables II and III. 

(r-jr 5- fo -xt 

L -I 

where X - arithmetic mean 

N - number of variates 



i-i 
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X - values of variate 


±- 


summation of 1 to N 


<r - standard deviation from the.' mean 
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TABLE II 


ARITHMETIC MEANS AND STANDARD DEVIATIONS 
(Not Used For Lag) 


East Components 


Variable 

Mean 

Standard Deviation 

Ae 

-1,03 

4.45 

X 1E 

-0.124 

3*18 

X 2E 

-0.25 

2.93 

^E 

-1.14 

3.62 

■%E 

-0.12 

3.32 

x se 

-0.70 

3.82 


-2.45 

8.42 

X7E 

-1.27 

6.75 

%E 

—1.35 

4.51 

^E 

0.00 

4.81 

X 10E 

-1.19 

4.29 

X 11E 

-0.06 

4.42 


North Components 


Variable 

Mean 

Standard Deviation 

% 

-0.02 

5.20 

*1N 

-0.42 

4.36 

X 2N 

1-1.38 

4.23 

X 3N 

-2.30 

4.19 

X 4N 

-0,14 

4.25 

•^N 

-0.25 

4.06 

^N 

-0.01 

7.99 

X 7N 

♦0.31 

7.23 

^SN 

♦1.24 

5.54 

^9N 

-0.42 

5.93 

X 10N 

-0.72 

5.43 

X 11N 

♦1.18 

5.23 
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TABLE III 


ARIHLKETIC MEANS AND STANDARD DEVIATIONS 
(Used For Lag) 


East Components 


North Components 


Variable 

Mean 

Standard Deviation 

Ae 

—1.03 

4.45 

?ie 

*0.47 

3.16 

Y 2E 

*0,07 

3.08 

*3E 

-1.33 

3.43 

Y 4E 

-0.76 

3.47 

Y 5E 

-0,05 

3.54 

Y 6E 

—1,55 

7.90 


-1.39 

6.26 

Y 8E 

-1.27 

4.76 

*9E 

-0,18 

4.88 

y ioe 

•0.83 

3.81 

y iie 

—0.59 

4.67 


Variable 

Lean 

Standard Deviation 


—0.02 

5,20 

Y 1N 

*0,22 

3,79 

Y 2N 

*0.61 

4.21 

Y 3N 

*0.20 

3.57 

Y 4N 

-0.75 

4.24 

Y 5N 

-0,28 

3,90 

Y 6 N 

♦0.57 

7.68 

Y 7N 

-0.24 

7.33 

Y 8N 

♦0.28 

5.50 

Y 9N 

-0.55 

5.74 

y ion 

-0.23 

5.25 

Y m 

-0.15 

5.27 
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TABLE IV 


TABLE OF SIMPLE CORRELATION COEFFICIENTS OF 
EASTWARD SPECIFIC MOMENTUM CHANGES AT 500 
MILLIBAR LEVEL TOTH EASTWARD ACCELERATION 
OF SURFACE UOYI 


Section 


Simple Correlation Coefficients 
No Lag Surface Low Lagging 1 Day 


South 

-.059 

+.082 

East 

4-. 203 

+.210 

North 

-.021 

+.038 

West 

-.013 

+.152 

Center 

+.206 

+.122 

Total All Sections 

+.144 

+.208 

Total - North,West,South,East 

+.094 

+.193 

North + East 

+.190 

+.114 

Yfest +■ South 

-.148 

+.209 

North + South 

+.041 

+.008 

West +■ East 

+.033 

+.301 
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TABLE V 


TABLE OF SIMPLE CORRELATION COEFFICIENTS OF 
NORTHWARD SPECIFIC MOMENTUM CHANGES AT 500 
MILLIBAR LEVEL WITH NORTHWARD ACCELERATION 
OF SURFACE LOW 


Section 


Simple Correlation Coefficients 


No Lag 

Surface Low Lagging 1 Day 

South 

<-.189 

-.120 

East 

-.141 

+.265 

North 

-.129 

‘ +.293 

West 

+.133 

+.018 

Center 

+.250 

-.081 

Total All Sections 

+.192 

+.201 

Total - North,West,South,East 

-.002 

+.224 

North + East 

-.204 

+.378 

West + South 

+.329 

-.077 

North + South 

+.109 

*,086 

West + East 

-.028 

+.225 
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Single correlation coeficients were then computed, using the 
accelerations of the surface low to the specific momentum changes at 
the 500 millibar level. In all 65 values of each variate were used. 
The correlation coeficients were computed from the formula 


(3) 




-- -i J: (xt-x) w 

(Tjc (7"y 


-Y) 


where 

x and y 
x and y 
Ox V 


s correlation coeficient 
s variates correlated 

s arithmentic means of variates correlated 
= standard deviations of the variates correlated. 


From (2) and (3) it can be seen that in order for the existence of a 
high correlation coefficient the respective signs of (X~X ) and (y~V ) 
must not only be consistent but their magnitudes must also be very 
nearly the same. The simple correlation coefficients are presented in 
Tables IV and V. Table IV gives correlations of eastward components and 

t 

Table V northward components. 

Next certain of the variates which yielded the highest values of 
simple correlations were selected to obtain multiple correlations. The 
multiple correlation coefficient is written 1.23) and read the 
multiple correlation of variate 1 on variates 2 and 3. The formula 
for the multiple correlation coefficient is 
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where 


•^OO 

Slsfrl) 

-Asfi/J y l'( }> ) ^(33) 

and R 11 is the minor of yo(it) • In (4) , the sO s are simple corre¬ 
lations between the variates. In general in order to get high multiple 
correlation the simple correlations must be high. 

The standard error of estimate for each multiple correlation was 
computed. The formula is, 

(5) / a r 

5 ^ C|rT (7 1 ) . 

It is called the standard error of estimating X^ from Xg and X^. 
Examining ( 4 ) and (5), it is seen that high values of multiple corre¬ 
lation will give low values of standard error of estimate. The value 
of the standard error of estimate is also directly proportional to the 
standard deviation of the dependent variate correlated. It is seen 
that the smaller the standard error of estimate the better the chance 
of predicting the dependent variate. The multiple correlations con¬ 
sidered and the corresponding standard errors of estimate are pre¬ 
sented in Table VI, 

Although the correlations are not high as correlation coefficients 
are generally rated in most fields, they do show a rather interesting 
pattern. 
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TABLE VI 


TABLE OF MULTIPLE CORRELATION COEFFICIENTS 
AND STANDARD ERRORS OF ESTIMATE 


Variates Correlated 

Correlation 

Coefficient 

’/>£ y»z) 

.393 

iab 'Xsb y„ £ 

.373 

(Ae 'Y/ib Yj.b) 

.373 

• XseX*e) 

.279 

^ (Ab ’Xsb Yie) 

.281 

^(Ab ‘X>e y^j 

.314 

^(Aa> * x*v M 

.460 

* Ysb y*») 

to 

CA 

• 

^(A» 'XimXsm) 

.A25 

Am 1 y*») 

.398 

-^(An' x?m yw 

.388 


Standard 

Error of Estimate 
4.09 
4.15 
4.12 
4.27 
4.27 
4.22 
4.62 
4.67 
4.71 
4.77 
4.79 
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The north sector when correlated concurrently with the northward 
acceleration of the surface low shows a negative correlation of -. 129 , 
but when the specific momentum changes are lagged by 1 day, the corre¬ 
lation coefficient becomes a f .293. Thus an increase of northward 
momentum in the north sector in the latter case could be interpreted 
as a surge of mass northward away from the position of the surface 
center. This surge of mass northward leaves deficit of mass to the 
north of the surface center and the effect is reflected the following 
day by an increase in the northward movement of the surface low. 

Consider now the west + south section’s values. The correlation 
of the northward component of this combination of sections to the 
northward acceleration of the surface low concurrently gives a value 
of . 329 , but when the south +■ west sections are lagged by one day 
the correlation coefficient becomes - .077. Thus a northward surge 
in the rear sectors causes an excess of mass to the south of the low 
and the effect is concurrently noted by a northward acceleration of 
the surface center. 

The array of correlations for the eastward components indicates 
that such an explanation as given above does not suffice for eastward 
accelerations of the surface low. The best simple correlations are 
obtained by the combination west +• east sectors lagged by 1 day and 
the east sector lagged by 1 day, the correlation coefficients being 
respectively +> .301 and r .210. Although the combination of west * 
east gives the greater correlation, its standard deviation is 4*67 as 
compared to 3.08 for the east section alone,which reduces its apparent 
advantage. 

The center sector which is analogous to the 500 millibar wind 
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shows significant correlation, only vdth no lag. The correlation 
coefficients obtained are smaller than those obtained by other in¬ 
vestigators correlating upper wind to movement of surface lows, but 
it must be remeribered that the correlations presented here are corre¬ 
lations of accelerations rather than velocities and hence it was ex¬ 
pected that they would be smaller. 

When the specific momentum changes at the 500 millibar lavel are 
lagged by the accelerations of the surface low the situation lends it¬ 
self to a forecasting use. Two regression equations are presented 
below. The regression equation is an implement for predicting the 
value of the dependent variate when the values of the independent 
variates are put into the equation. The general form of the regression 
equation for 3 variates is 

( 6 ) 

fin ^ O 

(77 <r> cn 

The symbols in equation (6) have been described in the preceeding pages. 
Regression equations for predicting both eastward and northward accel¬ 
eration of the surface low have been computed. The equations are 

(7) r 

As •-.A5’ y»M ~ h 7 / 2 - } JSAe) - ± V./Z- 

A m = /.3/ Yw +,3a2~y&/f i- } 3(4*) 

In (7) the independent variates are in units of b LoX./cla,y %. . 

T% * 
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Their corresponding multiple correlations of *• .373 and +- .438 
respectively are not high enough to give good quantitative results, 
although the qualitative indications should be significant. 

For qualitative prediction of northward acceleration of surface 
centers the north component of specific momentum change in the north 
section lagged one day was observed to be, with few exceptions, of the 
same sign (compared to respective means) as the northward acceleration 
of the surface low. 
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CHAPTER 17 


CONCLUSIONS 

It has been shown by this investigation that there is correlation 

% 

between accelerations of surface lows and specific momentum changes at 
the 500 millibar level. This, in effect, is also a correlation between 
day-to-day zonal and meridional index changes and movements of surface 
lows. It shows the effect of index on a micro-scale. 

Although the correlation coefficients obtained were comparatively 
low, this was not too surprising. There are several factors contributing 
to this: 

1. Only one level of the whole atmosphere is considered. If the 
investigation were to be conducted for numerous levels the results should 
be more significant. 

2. The correlations of accelerations rather than velocities leads 

to smaller numerical coefficients, especially if the day to day velocities 
of the surface low are fairly uniform. 

3. The distance of 600 nautical miles between the center and north, 
south, east, and west sections was arbitrarily chosen, and that it is 
the best distance to use is open to question. In several instances it 
was observed that specific momentum values were obtained from outside 
the limits of the 500 millibar trough reflected from the surface low. 

In cases where the east or west sections from one day to the next passes 
from the confines of the trough across the ridge line the changes in the 
north component of the specific momentum would be quite large and not 
necessarily reflected in accelerations of the surface low. It is sur¬ 
mised that using a shorter distance from center to outlying sections 
would give better results. 
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4* Geostrophic winds were used throughout. In regions of marked 
curvature the errors are quite significant. 

5. The scale of the maps (5) was quite small. Even though dividers 
were used, for wind velocities greater than 75 miles per hour errors of 
less than 10% are difficult to avoid. 

Although the correlation coefficients obtained in this investigation 
are smaller than those obtained in other investigations which correlated 
velocities, it does not mean that they are less useful. For when predict¬ 
ing from velocity correlations one must predict the total displacement, 
but in predicting from acceleration correlations only the expected change 
from the previous days displacement is predicted. In cases where the 
differences in day^to-day displacements of the surface low are small 
compared to the displacements themselves the value of predicting from 
acceleration correlations is enhanced. Thus in such cases the regression 
equations based on multiple correlation of acceleration and analogous 
specific momentum changes present a refinement over regression equations 
based on velocity correlations. 

This investigation has scarcely made a start in the possible lines 
of investigation. The best distances of center to outer section is not 
determined. The investigation could be carried out at other levels if 
data is available. The time of lag that gives best results is un¬ 
determined. The lag enforced by the data in this investigation was 32g 
hours. It is most probable that a different time lag would give better 
results. 

In conclusion, it is felt that the results obtained are significant 
and that micro scale changes in meridional and zonal indices are re- 
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fleeted in the day to day accelerations of surface lows. However, 
much more extensive work needs to be done to complete the problem. 


23 



BIBIDGRAPHY 


1. NAMAIS, Jerome, and CLAPP, Philip F., Studies of the Motion and 
Development of Long Waves in the Westerlies. Journal of Meteorology, 
V. 1, Nos. 3-4, PP. 57-77, December 1944. 

2. ROSSBI, C.- G., Relation between Variations in the Intensity of 
the Zonal Circulation of the Atmosphere and the Displacements of 
Semi-Permanent Centers of Action, J. Marine Res,, vol. 2, No. 1, 

PP. 38-54, 1939. 

3. ROSSBY, C.- G., "Climate and Man", Yearbook of Agriculture for 
1941, U. S. Government Printing Office, Washington, D. C. 

4. HAURWITZ, Bernhard, Dynamic Meteorology, PP. 265-267, Me Graw- 
Hill Book Company, 1941. 

5. Historical Weather Map Series, 500 millibar and Surface, U. S. Army 
Air Forces. 

6. WILLETT, H. C,, Final Report of the Weather Bureau - M.I.T. Ex¬ 
tended Forecasting Project for the Fiscal Year July 1, 1946 - July 
1, 1947. 

7. AUSTIN, J. H., An Empirical Study of Certain Rules for Forecasting 
the Movement and Intensity of Cyclones. Journal of Meteorology, 
4 : 16 - 18 . 

8. LONGIEY, B. W., A Study of the Relationship Between the 700 
Millibar Flow and the Movement of Surface Pressure Centers. 

Journal of Meteorology, 4:202-204. 

9. RUNK, T. W., A Study of the Relationship of Displacement of 
Cyclones to Various Aspects of the Upper Air Flow, U. S. Naval 
Postgraduate School, Annapolis, Maryland, 1948. 


24 



